EPJ manuscript No. 

(will be inserted by the editor) 



Non-Perturbative Contribution to the Thrust Distribution in 
e+e~ Annihilation 

R.A. Davison'' and B.R. Webber 

Cavendish Laboratory, J.J. Thomson Avenue, Cambridge CBS OHE, UK 
19 September 2008 

Abstract. We re-evaluate the non-perturbative contribution to the thrust distribution in e~ hadrons, 
in the light of the latest experimental data and the recent NNLO perturbative calculation of this quantity. 
By extending the calculation to NNLO+NLL accuracy, we perform the most detailed study to date of the 
effects of non-perturbative physics on this observable. In particular, we investigate how well a model based 
on a low-scale QCD effective coupling can account for such effects. We find that the difference between 
the improved perturbative distribution and the experimental data is consistent with a l/Q-dependent 
non-perturbative shift in the distribution, as predicted by the effective coupling model. Best fit values of 
as (91.2 GeV) = 0.11641E;;°E;^^ and ao (2 GeV) = 0.59 ± 0.03 are obtained with xVd-o.f. = 1.09. This is 
consistent with NLO-I-NLL results but the quality of fit is improved. The agreement in ao is non-trivial 
because a part of the 1/Q-dependent contribution (the infrared renormalon) is included in the NNLO 
perturbative correction. 

PACS. 13. 66. Be Hadron production in e~ interactions - 12.38.Cy Summation of QCD perturbation 
theory - 12.38.Lg Other nonperturbative QCD calculations 

1 Introduction 

One of the most common and successful ways of test- 
ing QCD has been by investigating the distribution of 
event shapes in e^e~ hadrons, which have been mea- 
sured accurately over a range of centre-of-mass energies 
(14 GeV < Q < 207 GeV), and provide a useful way of 
evaluating the strong coupling constant as. 

The main obstruction to obtaining an accurate value 
of as from these distributions is not due to a lack of pre- 
cise data but to dominant errors in the theoretical calcu- 
lation of the distributions. In particular, there are non- 
perturbative effects that cannot yet be calculated from 
first principles but cause power-suppressed corrections that 
can be significant at experimentally accessible energy scales 
In the case of the thrust distribution da/dT, previous work 
has shown that matching as with a low-scale effective cou- 
pling Q!eff which extrapolates below some infra-red match- 
ing scale Hi results in a 1 / Q-dependent shift in the distri- 
bution that accounts well for the discrepancy between the 
experimental and perturbative results W. 

The presence of l/Q corrections in event shapes is a 
generic expectation based on the renormalon analysis of 
perturbation theory, which implies an ambiguity of that 
order in the perturbative predictions for these observables 
(see [2J[3] for reviews) . The low-scale effective coupling hy- 
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pothesis [1] leads to universality relations between the cor- 
rections to different observables, valid to lowest order in 
the effective coupling, and to a well-defined prescription 
for matching the perturbative and non-perturbative con- 
tributions. 

The calculation of Ref. [T] was performed to NLO+NLL 
accuracy, i.e. terms up to O (a^) were retained exactly 
while exponentiating logarithmically-enhanced terms of 
the form ln"+\l - T) and a'; ln"(l - T) were summed 
to all orders. In the present paper, the recent evaluation of 
the NNLO term (i.e O (tts)) in the fixed-order perturba- 
tion series expansion of the thrust distribution [SJ[5] is used 
to refine the perturbative calculation of the distribution to 
NNLO-I-NLL accuracy and thus to reduce the uncertainty 
present in the theoretical prediction. A low-scale effective 
coupling is then introduced and matched to NNLO . This is 
again found to be a good method for dealing with the non- 
perturbative shift. By comparing the NNLO-|-NLL-|-shift 
results with the latest experimental distributions, values 
of as and 

1 f^" 

Qfo = — / d^i aoff (m) (1) 
Jo 

are obtained. These are consistent with those determined 
to NLO-I-NLL accuracy. The agreement is non-trivial be- 
cause a part of the 1 /Q-dependent contribution - the in- 
frared renormalon - is included in the NNLO perturbative 
correction. 

The organisation of the paper is as follows. In Sect. [2] 
we briefiy recall the relevant properties of the thrust distri- 
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bution, the fixed-order calculation and the resumniation of 
large logarithms. Sect.[3]presents the predictions of pertur- 
bative NNLO+NLL matching and the power dependence 
of the discrepancy with experimental data. The matching 
to the low-scale efi^ective coupling and comparisons with 
data are performed in Sect. |4j and our conclusions are 
presented in Sect. [5j 

2 Perturbative calculation of the thrust 
distribution 

We recall that the thrust T is a measure of the distribution 
of momenta of the final state hadrons: 



T ■ 



E^i \pl-n\ 



(2) 



where "n is a unit vector and we sum over the 3-momentum 
of each final-state hadron in the centre-of-mass frame. 
Theoretical calculations of thrust are performed by sum- 
ming over the individual final state partons, as the hadro- 
nisation process is still not well understood. T can vary 
between the limits T = 1 for back-to-back jets and ^ = | 
for a uniform angular distribution of hadrons. 

For comparison with experiments, it is the thrust dis- 
tribution 

which is relevant, where a is the total cross-section for 
e^e~ — !■ hadrons. In calculations it is more convenient to 
use the event shape variable 



t EE 1 - T, 



(4) 



which has the two-jet limit t — Q. The distribution away 
from this limit therefore depends directly upon the pro- 
duction of extra final-state partons at QCD vertices, and 
hence is ideal for testing QCD and evaluating a^. The 
normalised thrust cross section is then defined as 



Jo a dt Ji_t a 



1 da 
dT' 



(5) 



2.1 Fixed-order calculations 



The perturbative expansion of the normalised thrust cross 
section has the general form 



R{t) = 1 -1^ OsRiit) + alR2{t) + alRz{t) + 



(6) 



where i?i {t) is the leading order (LO) coefficient, i?2 {t) is 
the next-to- leading order (NLO) coefficient, i?3 {t) is the 
next-to-next-to-leading order (NNLO) coefficient etc. and 
OLs = as/2TT. Solving the renormalisation group equation 
for the running coupling to NNLO gives 



as (Mi?) 



2tt 



1 



/3i In L 1 



InL 



Po 



(7) 



where fiji is some chosen renormalisation scale (we take 
fJ-B. = Q except where stated otherwise). 



Po 

P2 



WN ~2Nf 



432 



6 

(28577V^ 



- ^NNf - iGpNF 



MCINf 



6 

79 NNl) , 



(8) 

with Cf {N^ - 1)/2N for an SU(iV) gauge theory with 
Nf active ffavours {N = 3 for QCD and Nf = 5 at all 



energies considered here) and L = ln(/i'^/yl^^'^), --^^^ 
being the 5-ffavour QCD scale in the modified minimal 
subtraction renormalisation scheme. 

A numerical Monte Carlo program, EERAD3 [7], has re- 
cently been developed which computes the process e~^e~ 
jets to NNLO in as via the decay of a virtual neutral gauge 
boson (7 or Z") to between three and five partons [5J[5]F] 
The EERAD3 predictions for the thrust distribution at a 
variety of centre-of-mass energies Q spanning the range 
14 GeV to 206 GeV are shown by the green/lighter curves 
in Figs. T]|3 The values of as (Q) were calculated using 



A. 



(5) 



i{5}_ 
MS 

as (91.2 GeV) 



= 0.204 GeV, corresponding to the world average 



0.1176 ig. 



2.2 Resummation of large logarithms 



The enhancement of the distribution at low t due to soft 
or coUinear gluon emission (as seen in Figs.[l]-[3| is present 
at all orders in perturbation theory: the dominant term at 
nth order is typically of the form 



1 da 
a dt 



1 



■In 



2n-l 



(9) 



Thus we see that at low t the condition as ^ 1 is not 
sufficient for a fixed-order prediction in perturbation the- 
ory to be accurate. Instead, we require UgL^ <C 1, where 
L = ln(l/t). To obtain accurate predictions in the two- 
jet limit t — > 0, we must therefore take account of these 
enhanced terms at all orders in perturbation theory by 
resumming them. 

Resummation of large logarithms is possible for event 
shape variables y that exponentiate [10,, i.e. their corre- 
sponding normalised cross section can be written in the 
form 



R{v)^C (as) S (y, as) + D (y, as) , 



(10) 



^ A recent calculation |8j finds some discrepancies with 
Refs. 5,6 , but these are not significant in the kinematic re- 
gions that we consider. 
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NNLO 

NNLO+NLL 

TASSO(14GeV) i — i — 

Q=14 GeV 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 



NNLO 

NNLO+NLL 

TASSO (35 GeV| i — i — 
JADE (35 GeV) i— ^ 
0=35 GeV 
"3=0.1375 




NNLO 

NNLO+NLL 

AMY (54.5 GeV) i — i — 
L3 (55.3 GeV) i— ^ 

Q=55 GeV 
"3=0.1273 




NNLO 

NNLO+NLL 

L3 (65.4 GeV) i — i — 
DELPHI (66.0 GeV) i— ^ 

Q=66 GeV 
"3=0.1236 







Fig. 1. Fixed-order (NNLO), resummed (NNLO+NLL) and 
experimental thrust distributions: Q — 14 — 66 GeV. 



Fig. 2. Fixed-order (NNLO), resummed (NNLO-fNLL) and 
experimental thrust distributions: Q — 91 — 183 GeV. 
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NNLO 

NNLO+NLL 

ALEPH jiag^O GeV) i * i 
DELPHI (189.0 GeV) i x ' 

Q=189 GeV 
"s=0.1060 




NNLO 

NNLO+NLL 

L3(200 GeV| i — i — i 
ALEPH (200 GeV) i * i 
DELPHI (200 GeV) i x i 
0=200 GeV 
"s=0.1052 




Fig. 3. Fixed-order (NNLO), resummed (NNLO+NLL) and 
experimental thrust distributions: Q = 189 — 207 GeV. 



where 



C(a,) - 1 + 



In 17 (j/, as) 



oo n-f-l 

EE 

n— 1 m—1 



G-n r 71 



Lgi (a^L) + g2 {oisL) + a^g^ {a^L) + 



L — ln(l/?/) and D(y,as) is a remainder function that 
vanishes order-by-order in perturbation theory in the two- 
jet limit y 0. The functions gi {agL) are power se- 
ries in UgL (with no leading constant term) and hence 
Lgi (agL) sums all leading logarithms a"L"+^, g2{asL) 
sums all ncxt-to-leading logarithms (NLL) a"L" and the 
subdominant logarithmic terms d^U^ with < to < n 
are contained in the gs, g^, . . . terms. The functions gi thus 
resum the logarithmic contributions at all orders in per- 
turbation theory, and knowledge of their form allows us 
to make accurate perturbative predictions in the range 



agL < 1 - a significant improvement on the fixed-order 
range a^L^ <C 1. 

For thrust, the first two functions can be determined 
analytically by using the coherent branching formalism 
I12| . which uses consecutive branchings from an initial quark- 
antiquark state to produce multi-parton final states to 
NLL accuracy. The results of this calculation depend upon 
the jet mass distribution j[Q^,k^^ - the probability of 
producing a final state jet with invariant mass from a 
parent parton produced in a hard process at scale - 
and its Laplace transform (Q^) • To the required accu- 
racy, the thrust distribution is 



1 da 
a dt 



91 

2TTi 



(12) 



where the contour C runs parallel to the imaginary axis 
on the right of all singularities of the integrand. 



inJ/;(g2 



1 du 



100 



l-K 



K = N 



as (m)\~^ , 



"'^'d/iV a,(/i) 



27r 



67 
18 



TT 

"6 



-Nf 



(13) 



(14) 



This expression demonstrates explicitly that the diver- 
gence of tts (/i) at low ^ will affect the perturbative thrust 
distribution - such effects are related to the renormalon 
mentioned earlier. To NLL accuracy, however, we can ne- 
glect the low ^ region (although we will return to it in 
Sect. |4]) to give the thrust resummation functions 



gi {usL) = 2/i {PqUsL) , 

32 {asL) = 2/2 (/3oa,i) - \nr[l - 2/i (/3oa,L) 
-2/3oasL/a/?oa.i)], 



(15) 



where 



[(l-22;)ln(l-2x) 



-2(1 - .t) ln(l - a;)] 
CfK 



PI 
3Cf 



[21n(l 



ln(l -x) 



-ln(l 
2Cf7e 



-2x)] 
[ln(l-a;) 



-ln(l-2a;)]-^^[ln(l-2x) 



21n(l 



1 



ln{l-2x) -ln^(l 



(16) 



By writing the K dependence in the form shown in ( 13 1, we 



change from the MS renormalisation scheme to the so-called 
bremsstrahlung scheme [13j . 
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with r the Euler r'-function, the Euler constant, and 
Cp, K and /3„ the constants previously defined. 

By combining these with the fixed-order calculation, 
we can obtain a new estimate of the normalised cross sec- 
tion to NLL accuracy. This should particularly improve 
the fixed-order estimate in the two-jet region, where L 
becomes large. Naively we would simply calculate R (t) as 
defined in Eq. (10), but it turns out to be considerably 
simpler to consider lni?(t), as we recall next. 



2.3 Log-R matching 



In the log-R matching scheme, we rewrite the exponenti- 
ation formula as 

In R{t)^ F (as) +lnE {t, as) +H{t,as), (17) 

where F (as) is a power series in and H (t, ag) denotes 
the remainder function which vanishes as t 0. 

For a fixed-order perturbative calculation of R (t) to 
order M, we can write Eq. ^ as 



In R (t) = In ^1 + a'^R,, (t)^ 

M 1 / \ ^ 

= ^«';i?„(i)--K]«:i?„w (18) 

n=l \n=l ) 



The matched estimate is obtained by combining the Mth 
order perturbative result with the resummed contribu- 
tions and subtracting the terms of order < M in In J7 (as 
these are already accounted for in the fixed-order terms). 
Thus for a fixed-order calculation to order a^, the matched 
estimate after resumming large logarithms to NLL accu- 
racy is 

In R (t) = Lgi {usL) + 52 {asL) 

+ as (i?i {t)-GnL-Gi2L^) 

+ al (t) - I [Ri it)f - - G23L' 

+ al (^i?3 (t) - Ri (t) i?2 (t) + I [Ri {t)f 

(19) 

The coefficients Gnm can be extracted by expanding 
the functions gi {asL) and 52 (ctgL) as power series in asL 



and comparing them with the definition (111 of G 



Gil — SCp, 
G12 — — 2Cf, 



G22 
G23 
G33 



G 

UStt^Gf + (169 - 127r2) N - 22Nf] 
36 '■ ^ ' ■' 

Cf 



{UN -2Nf), 
[2304C (3) C| - 7927:'^ NGf 



G 



34 



108 

(3197 - 1327r2) + (IO8 + lUn^) GfNf 
(1024 - 247r2) NNf - 68A^|] , 

-^GF{nN-2NF)\ 



(20) 



where C (3) = 1.202057.... 

There are two reasons why it is simpler to use this 
log-i? matching scheme rather than R matching (i.e. eval- 
uating Eq. (10) explicitly to NLL precision). Firstly, we 
do not have tobe concerned with the C (ag) and D {t, ag) 
terms in (10), for which we do not have analytic expres- 



sions but which contribute to the fixed-order calculation 
- these are contained in Ri{t), R2{t), etc. Secondly, it 
is easier to impose physical boundary conditions on the 
normalised cross section, namely 



R (^t — tjjifix) — 1: 

by definition of the normalised cross section, and 



dR 
~dt 



(21) 



(22) 



as there is an upper kinematic limit tamx on the thrust 
for a given number of final-state partons. Although the 
resummed logarithmic terms are small at high t, dR/dt 
is also small and so these terms can cause relatively large 
unphysical effects if we do not impose these conditions. 

The above constraints are automatically obeyed by 
the fixed-order terms i?„ {t) but not by the resummed 
terms, as we have neglected the subdominant logarithms 
53 (ttsi), g/i{asL) etc. To satisfy these constraints, we 
therefore require 

Q (t) = Lgi (asL) + 52 i^sL) ~ as {GuL + Gi2L^) 

- al {G22L^ + G23L^) - al (GgaL^ + 034^^) 

(23) 

and its first derivative to vanish at i = tmax- Q (t) cor- 
responds to the resummed logarithmic terms of order 
and higher and hence at small L, 



dt 



By making the replacement 



L ^ L = In 1 



(24) 



(25) 
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the boundary conditions are satisfied as L (tmax) = 0. This 
does introduce corrections to the expression for In R (t) 
but these are power-suppressed at small t: 



L (t) In 



In 1 - 



t 



+ t 



L{t)+[t 



and so L (t) L (t) in the important limit t ^ 0. 



(26) 



3 Results of NNLO+NLL matching 

To perform the matching, the integrated perturbation se- 
ries coefficients are required as in Eq. (19 1. For Ri (t), the 
analytic result is 



Ri (t) 



where 



■In" 



3 
10 

y 



t 



i-t 

8t + 6t^ 



Li2 (z) = 



4(1- 2t) In 



t 



1 - 2t 



16 



Li2 



1-t 



° , In (1 - a;) 
ax 



47r^ 
9 



(27) 
(28) 



is the dilogarithm function. i?2 (t) and R3 {t) were ob- 
tained by interpolating the differential results from EERAD3 
and then numerically integrating them. For R3 (t), the 
EERAD3 results were first smoothed by taking 



dR3 
dt 



iU) 



dR3 
dt 



dRa 
dt 



(29) 

repeatedly until a smooth curve was obtained. The peak 
near t = had to be reintroduced by hand, as this smooth- 
ing technique always results in the peak value being re- 
duced. 

R (t) was computed to NNLO-I-NLL precision using 
Eqs. (19 1 and (26) with tmax = 0.42 in L, as this is the 
maximum value of t kinematically allowed in the five par- 
ton limit. The differential cross section was then obtained 
by numerically differentiating R{t) . The results at a range 
of energies are shown by the red/darker curves in Figsll]l3j 
The values of as {Q) were calculated as described earlier 
for the unresummed NNLO (green/lighter) curves. The 
shaded area around each line shows the renormalisation 
scale uncertainty found by taking fij^ G [Q^/2, 2(3^] . 



Experiment 


Q/GeV 


Ref. 


No. Pts. 


2 

X 


lAbbO 


14.0 


[HI 


4 


8.2 


lAbbO 


22.0 


en 


6 


2.8 


TASSO 


35.0 




8 


0.7 


JADE 


35.0 


mi 


10 


10.5 


T 


41.4 




8 


3.4 


JADE 


44.0 


nn 


10 


3.8 


lAbbO 


44.0 


mi 


8 


6.8 


DELPHI 


45.0 


im 


11 


11.6 


AMY 


54.5 


M 


4 


4.9 


L3 


55.3 


M 


8 


3.2 


L3 


65.4 


ESI 


8 


7.5 


DELPHI 


66.0 


[13 


11 


14.5 


L3 


75.7 


IE] 


8 


1.9 


DELPHI 


76.0 


[13 


11 


10.3 


L3 


82.3 


M 


8 


4.0 


L3 


85.1 


M 


8 


3.6 


OPAL 


91.0 


M 


5 


11.9 


ALEPH 


91.2 


M 


27 


16.1 


DELPHI 


91.2 


[13 


11 


18.8 


SLD 


91.2 


M 


6 


2.7 


L3 


130.1 


[m 


10 


14.6 


ALEPH 


133.0 


M 


6 


7.2 


OPAL 


133.0 


M 


5 


6.5 


L3 


136.1 


[SI 


10 


37.3 


ALEPH 


161.0 


M 


6 


5.5 


L3 


161.3 


M 


10 


4.0 


ALEPH 


172.0 


M 


6 


14.0 


L3 


172.3 


[SI 


10 


2.1 


OPAL 


177.0 


M 


5 


1.1 


L3 


182.8 




10 


2.7 


ALEPH 


183.0 


M 


6 


4.0 


DELPHI 


183.0 


[13 


13 


33.1 


L3 


188.6 


[m 


10 


3.4 


ALEPH 


189.0 


M 


6 


6.7 


DELPHI 


189.0 


[13 


13 


22.7 


DELPHI 


192.0 


[13 


13 


12.1 


L3 


194.4 


M 


10 


1.2 


T~\T7'T TITTT 

DELPHI 


196.0 


[13 


13 


39.7 


r\Ti A T 

U-tAIj 




1-1 ol 


r 



inn 


ALEPH 


200.0 




6 


21.0 


DELPHI 


200.0 


[13 


13 


7.1 


L3 


200.0 


M 


9 


6.5 


DELPHI 


202.0 


[13 


13 


14.9 


DELPHI 


205.0 


[13 


13 


12.6 


ALEPH 


206.0 


M 


6 


7.0 


L3 


206.2 


M 


10 


10.0 


DELPHI 


207.0 


m 


13 


11.7 


Total 






430 


466.0 



Table 1. Data sets used and best-fit contributions. 



3.1 Comparison with experimental data 

The matched, resummed differential thrust distribution 
was compared with data from a wide range of experiments, 
as listed in Table [l] The points in Figs. l]|3 show the data 
at an illustrative selection of energies. The error bars rep- 
resent the experimental statistical and systematic errors, 
added in quadrature. 



There are a few features common to the graphs at 
all energies. Firstly, the resummed distribution and the 
NNLO distribution are almost identical away from the 
two-jet region. However, in this low-< limit the resummed 
distribution peaks, in line with the experimental data, 
whereas the NNLO distribution carries on increasing. Thus 
resummation has significantly improved the theoretical 
prediction in the two-jet limit, as we had expected. 
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It should be noted that the kink around t — 0.33 in 
all of the curves is due to the LO term vanishing here 
for kinematic reasons. One would expect that with many 
higher-order perturbation theory terms taken into account 
(i.e. more partons present in the final state), this would 
gradually smoothen, in line with the experimental data. 

At all energies, the overall shape of the theoretical dis- 
tribution is similar to that of the data, but is shifted to 
a lower value of t. This apparent shift St has a clear en- 
ergy dependence - at the upper end of the energy range 
considered here, the NNLO-I-NLL and experimental dis- 
tributions are fairly close and the shift 5t is a very small 
correction. On decreasing the energy, the shift becomes 
more pronounced and at low energies the theoretical dis- 
tribution is clearly not consistent with the data. There is 
no obvious way that this could be remedied by the inclu- 
sion of sub-leading logarithms or higher fixed-order terms, 
and so we now turn to considering non-perturbative effects 
for an explanation. The increasing discrepancy at low en- 
ergies is also consistent with this interpretation, as we ex- 
pect such effects to have a 1/Q dependence, as mentioned 
in Sect. [T] To verify that these discrepancies are due to 
non-perturbative effects, the exact form of their energy 
dependence was investigated. 



3.2 Power dependence of discrepancies 

As both the experimental data and EERAD3 results are 
given as histograms, and not in terms of individual val- 
ues of t, the integrated thrust distribution R (t) should be 
slightly more accurate than da/dt as it does not involve 
the assumption of a uniform distribution over the width 
of each histogram bin At. 

Graphs of In (i?thcory — ^cxpt) against In (Q/GeV) were 
plotted for 0.025 < t < 0.24 and, anticipating corrections 
proportional to an inverse power of Q, a linear fit was 
made to each plot such that the gradient n of the straight 
line gives the power dependence of the required correction 

g[6) This t range 



(oc Q"). The results are shown in Figs, 
was chosen since at lower values of t there is no obvious 
straight line (due to the distributions peaking), and at 
higher values of t the percentage errors on the gradient 
become large due to i?thcory — -Roxpt quickly decreasing to 
zero (as both normalised cross-sections converge to 1). 

Although not totally conclusive, these results are con- 
sistent with power corrections of the form 1/Q, and we 
turn now to considering the quantitative form of these 
non-perturbative corrections to the thrust distribution. 



4 Non-perturbative corrections 

4.1 The low-scale effective coupling 

Although there are various ways to phenomenologically 
treat non-perturbative effects in QCD, one of the most in- 
tuitive is by means of a low-scale effective coupling [T]. In 
this approach, the running coupling ^ is replaced by an 
effective coupling acff (m) , which differs from the standard 




t = 0.025 

n = -0.85 ±0.08 




t = 0.05 

n = -1.45±0.1 




= 0.06 
= -0.88 + 0. 



DELPHI ^- 
ALEPH K 
OPAL K 
Bestnt - 





Fig. 4. Power dependence of corrections required to resolve 
theory/data discrepancy: t — 0.025 — 0.08. 
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3.B 4 4.2 4.4 4.6 4.B 5 5.2 5.4 3.5 3.6 3.7 3.B 3.9 4 4.1 4.2 4.3 4.4 4.5 4.6 

InlQ/GeVI ln(Q/GeV) 




2.5 3 3.5 4 4.5 5 5.5 2.5 3 3.5 4 4.5 5 5.5 6 

InfQ/GeVl ln(Q/GeV) 




-5 - - -5 



DELPHI h 
ALEPH h 
JADE h 



DELPHI >- 
ALEPH >- 

JADE 
TASSO — - 




Fig. 6. Power dependence of corrections required to resolve 
theory/data discrepancy: t = 0.18 — 0.24. 



perturbative as(/i) in the infra-red region where the latter 
diverges. Using this finite effective coupling allows us to 
use the formalism of perturbation theory to describe non- 
perturbative effects which cannot be probed using stan- 
dard perturbative QCD. 

Various forms for aeff(/i) have been proposed [22ll23] 
that have high-energy behaviour consistent with but 
we will not be concerned with their details here. The only 
parameter we require is the 'average' value of the effective 
coupling below the infra-red matching scale /i/ where 
and aoff begin to differ: 



f f^^' 

ao im) = — / dfi Qfcff in) . 
fJ-i Jo 



(30) 



infQ/Gev) We make the additional assumption that aoff is small 

Fig. 5. Power dependence of corrections required to resolve enough in the infra-red region that we can neglect terms 
theory/data discrepancy: t = 0.09 - 0.16. of order a'^g and higher. 
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4.2 Non-perturbative shift in thrust distribution 

In deriving the form of the NNLO+NLL prediction used 
earlier, the low /x region in Eq. (13 1 was neglected as 



it produced a subleading contribution. We now include 
this region by subtracting the fixed-order NNLO contri- 
bution from /i < Ai/ and replacing it with a contribution 
due to the effective coupling. We are thus removing the 
renormalon contributions to the perturbation series (up to 
NNLO) and incorporating all l/Q-dependent behaviour 
into aeff- 



Firstly, we note that the order of integration in (13 1 
can be changed, to give 



^ du 



27r 



1^ u 
Q2 



(31) 



Inserting the NNLO perturbative running coupling 
as in) = as (Mi?.) + (a^a) — In — 

IT 



(ma) 



In ^ In — 



(32) 



expanding the exponential to first ordeij^and integrating 
over the range < /i < /i/ gives an NNLO contribution of 



2Cf 
"s i^^R) 



(Mi?) 



2 / \ f^O 
(Mi?) — 



In 



Mi? 
Mi 



K 

Wo 



K 

Wo 



In- 



2 Mi? 
Mi 



Wo 



In 



Mi? 
Mi 



(33) 

It should be noted that t is the conjugate variable to 



in the Laplace transform ( 12 1 and thus the first-order 



expansion of the exponential will only be a valid approxi- 
mation in the limit i ^ Mi/Q- Below this, we would need 
to retain higher order terms in the expansion, which would 
require us to have a specific form for ctcff (m)- 

Following a similar procedure with (m) in the place 
of as (m) gives a non-perturbative contribution of 



2Cf 



dfi acs (m) i^Q 



2Cf Mi 



ao(Mi)'^Q', (34) 



where we have neglected terms of order a'^g as previously 
noted. 

By adding this, after subtracting the perturbative con- 
tribution (33 1, we obtain the change in the quark jet mass 
distribution caused by changing from a perturbative to an 
effective coupling in the low-scale region below /i/. 



Higher-order terms in the expansion would give corrections 
of order which we neglect. 



Substituting the result into Eq. (12 1, we see that the 



effect of this non-perturbative contribution is to shift the 
thrust distribution by an amount St, such that 



1 da 
a dt 



lda_y'- 
a dt J 



(35) 



t+st 



where 



6t^- ^ < ao (Mi) - (Mi?) 

TT Q y 

2( N /^O A Mi? , ^ , 

- ots (Mi?) — In — + + 



2 r 



(Mi?) 




(36) 



to NNLO. This l/Q-dependent shift is precisely what is 
required to account for the differences between the per- 
turbative and experimental distributions seen in Sect. [3j 



4.3 Determination of otg and ckq 



By applying the shift (36) to the perturbative results, we 



expect to reduce significantly the differences between the 
theoretical and experimental distributions. Comparison of 
these differences to the predicted form of 6t allows us to 
estimate the values of a^ and as ■ 

Maximum accuracy was obtained by comparing the 
experimental distribution with a discretely-defined theo- 
retical distribution 



At) 
At 



R{t) 



(37) 



where At is the bin width of the experimental data. 

The NNLO-|-NLL-|-shift distribution was calculated as 

a function of an and L. This calculation was performed 

" MS 

for 0.05 < t < 0.33, at the centre-of-mass energies listed 
previously in Table [l] (i.e. in the range 14 GeV< Q < 207 
GeV). was calculated for each pair of input parameters, 
with its minimum corresponding to the best-fit values. 

The upper limit for the fits was chosen as t = 0.33 
since the difference between the theoretical and exper- 
imental distributions above this value is largely due to 
the small number of final state partons in the theoret- 
ical calculation, as previously explained, rather than to 
any non-perturbative effects. We noted previously that 
the non-perturbative results are strictly valid only in the 
range t 3> ^ii/Q', in fact we found satisfactory fits using 
an energy-dependent lower cut-off t > maxj/ij/Q, 0.05}. 

For infra-red matching scale jij ~2 GeV, best fit val- 
ues of 

ao (2 GeV) = 0.59 ± 0.03 , 

yi(5) - iq0+° °25 - ^^^^ 



GeV 
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A(51 



Fig. 



7. X contour plot in ( , qq ) space. 



were obtained, with x^/d.o.f. ^ 466.0/428 « 1.09. The 
quoted errors correspond to one standard deviation, com- 
puted as recommended by the Particle Data Group [9] : the 
value of corresponding to the la (68.3% C.L.) contour 
was rescaled by the value of x^/d.o.f., giving — 480.6, 
i.e. Ax^ = 14.6. 

The contribution to from each data set is shown in 
Table [T] It should be noted that the few data sets with 
X^/no. pts. 3> 1 are not generally inconsistent with the 
shifted distribution, but simply have a few outlying points 
giving a large contribution. 

The contour plot in Fig. [7] shows the ranges of ao and 
^Ws ^^^'"^ within Ax^ of the best-fit value of x^, 

and also demonstrates the correlation between these two 
parameters. 

Varying the renormalisation scale //^ £ 2Q^] 
gave best fit values in the range (2 GeV) = 0.585, 
^ 0.173 GeV to ao (2 GeV) = 0.598, = 0.210 

IVl O 111 ij 

GeV with no significant change in the quality of fit. Thus 
we find 



A. 



(5) _ 
MS 



0.190 



-0.025+0.020 
-0.022-0.017 



GeV 



(39) 



where the first error is the combined experimental statis- 
tical and systematic error and the second is due to the 
theoretical renormalisation scale uncertainty. The corre- 
sponding strong coupling constant is 



as (91.2 GeV) = 0.1164; 



-0.0022+0.0017 
-0.0021-0.0016 ' 



or, combining all the errors in quadrature. 



as (91.2 GeV) = 0.1164; 



-0.0028 
-0.0026 ' 



(40) 



(41) 



in good agreement with the world average value of 0.1176 [S] 
To assess the importance of the NNLO terms, the anal- 
ysis was repeated with all those terms omitted, i.e. com- 
bining NLO-I-NLL in perturbation theory with Eq. ( 36 1 



without the O^a";) contribution. The resulting best fit val- 
ues were 

ao (2 GeV) = 0.51 ± 0.04 , 



A 



(5) 
'MS 



0.214l°:^f,iro^2^ GeV , (42) 



as (91.2 GeV)-0.1185™^H 



-0.027-0.026 
+0.0025+0.0027 



with xVd.o.f. = 515.1/428 w 1.20. Thus the NLO and 
NNLO results are consistent but the inclusion of NNLO 
terms consistently in both the perturbative prediction and 
the power correction improves the quality of the fit and 
reduces the errors. 

The most complete previous NLO study along similar 
lines 24J, combining NLO-f-NLL in perturbation theory 
with the NLO equivalent of Eq. ( 36 ) and covering a variety 



of event shapes but a slightly narrower range of energies 
than that used here, obtained the overall best fit at 



as (91.2 GeV) 



n 1 1 71 +0.0032 
U.ii/i_Q^0020, 



ao (2 GeV) - 0.513^^011 



(43) 



in good agreement with our results. Their fit to the thrust 
distribution alone gave 



as (91.2 GeV) = 0.1173+° °°^^ 



-0.0051' 



ao (2 GeV) = 0.492t^;°«4 



(44) 



).070 



also in good agreement. 

In the recent NNLO analysis [3S], a range of event 
shapes at energies at and above 91.2 GeV were fitted 
without resummation; non-perturbative effects were es- 
timated using Monte Carlo event generators. The value 
obtained for the strong coupling was as (91.2 GeV) — 
0.1240 ±0.0033. 

To estimate the dependence of our results upon the 
infra-red matching scale, a fit with /i/ = 3 GeV was made, 

yielding ao (3 GeV) = 0.458±0.025 and = 0.202^°;°^^, 

with x^/d.o.f. w 1.09. Thus the fit remains good and the 

f 5) 

value obtained for A^^ is stable under variation of /i/, 
while the value of ao decreases as expected for a running 
effective coupling. Indeed, the implied mean value of a^s 
in the range 2-3 GeV, 

aoff = 3ao(3 GeV) - 2ao(2 GeV) = 0.19 ±0.10 (45) 

is consistent with the perturbative value as (2.5 GeV) = 
0.26. 



4.4 Final comparison with experimental distributions 



Figures [8|T0| show the final (NNLO+NLL+shift) theoret- 
iistributio 



ical distributions in comparison to the experimental ones, 
with the best-fit values of ao and as assumed. The shaded 
area around the unshifted distribution is the renormalisa- 
tion scale uncertainty found by varying fij^ e [Q^ /2, 2(5^] , 
and the shaded area around the shifted distribution is 
the corresponding error found by varying between the 
best fit limits obtained previously (ao (2 GeV) = 0.585, 



= 0.173 GeV and ao (2 GeV) 



0.598, yl^= 0.210 



0023 



GeV). 

It is clearly seen that inclusion of the shift results in a 
significantly more accurate distribution over the fit range, 
particularly for the lower energies. As the best fit value of 
as is very close to the world average, the unshifted distri- 
butions here are essentially the same as those in Figs. T][3 




Fig. 8. Comparison of shifted, unshifted and experimental Fig. 9. Comparison of shifted, unshifted and experimental 
thrust distributions: Q = 14 — 66 GeV. thrust distributions: (5 = 91 — 183 GeV. 
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Unshilted 

Shifted 

ALEPH jiag^O GeV) i— • — i 
DELPHI (189.0 GeV) i ic 
Q=189 GeV 
"S=0.1050 




Unshifted 

Shifted 

L3 (200 GeVI i — i — i 
ALEPH (200 GeV) . — i 
DELPHf (200 GeVI i ic 
Q=200 GeV 
"s=0.1042 




Fig. 10. Comparison of shifted, unshifted and experimental 
thrust distributions: Q = 189 - 207 GeV. 



The agreement of the ao and as values from the analy- 
sis at NNLO+NLL with those obtained at NLO+NLL is a 
non-trivial test of the low-scale effective coupling hypoth- 
esis. The presence of the 0{a^) term in Eq. (36 1, which 
amounts to about 80% of the 0{a1) term, means that 
we are not simply adding a l/Q correction to the pertur- 
bative result, but rather that we are regularizing the di- 
vergent renormalon contribution by modifying the strong 
coupling at low scales. This implies that the explicit non- 
perturbative 1 /Q shift applied to the perturbative predic- 
tion becomes smaller as higher orders are computed, and 
would eventually change sign at sufficiently high orders, 
as the renormalon contribution grows indefinitely. 

A similar analysis to that in this work could be re- 
peated for other event shape variables whose distributions 
have been determined perturbatively to NNLO and for 
which resummation of large logarithms is possible. Pertur- 
bative resummed calculations of such distributions have 
been performed [26 but non-perturbative effects have not 
been included in the way advocated here - they are not 
necessarily simple shifts as in the case of thrust. It would 
also be of interest to combine the present approach to 
non-perturbative effects with soft-coUinear effective the- 
ory, which permits the resummation of next-to-next-to- 
leading logarithms [27] . 
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5 Conclusions 

We have seen that the extension of the NNLO perturba- 
tive distribution to NNLO-I-NLL accuracy results in an 
improved matching with experiment, particularly in the 
low t region. 

Analysis of the difference between the perturbative 
and experimental distributions over a range of energies 
showed that l/Q power corrections were required to ac- 
count for this difference. Replacement of the perturba- 
tive strong coupling with an effective coupling below an 
infra-red matching scale was used to include such non- 
perturbative corrections in our theoretical calculation and 
resulted in a l/Q-dependent shift in the distribution. With 
best-fit values ao (2 GeV) = 0.59±0.03 and (91.2 GeV) = 
0.1164^q;qq26i this gave a significantly improved matching 
with the experimental distributions in the range 14 GeV < 
Q < 207 GeV. These values are consistent with those 
achieved in similar analyses to NLO, as well as with the 
world-average value of a^. 
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